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Abstract 

In this paper we investigate the macroscopic quantum tunnelling and the phase coherence prop- 
erty of the rotating Bose-Einstein condensates in both static and dynamic cases by using the mean 
field theory. 
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I. INTRODUCTION 



It is well known that the phase coherence and the macroscopic quantum phenomenon are 
among the important characteristics of the Bose-Einstein condensates ( BEC ), and have 
induced a lot of theoretical and experimental interests. Now, we came to know that the 
phase coherence property and the macroscopic tunnelling are so closely related that such 
a system was found to exhibit quite rich physics. On the other hand, motivated by the 
superfluidity of the dilute Bose gas, much work are focused on the rotating BEC recently 
||. In this paper, we investigate the macroscopic quantum tunnelling in the rotating BEC 
in various parameter ranges and in both static and dynamic cases with the help of the mean 
field approach, and the phase coherence property in this system will be emphasized. 

In Section II, a general description of our system and the Hamiltonian to begin with 
are given. In the next two sections, we focus on the static property when the rotating fre- 
quency is fixed. SU(2) coherent state path integral method is used in Section |T| and we find 
that if the parameters satisfy a specific condition, the ground state represents the super- 
position of two macroscopic quantum states. The energy splitting due to the macroscopic 
coherent tunnelling is also calculated by the instanton technique. In another parameter 
range, fluctuations of the phase and the particle number of the ground state are given in 
Section [IV] and it follows that the condition of the Josephson effect in this case is satisfied. 
In Section M, we consider the case where the rotating frequency varies with time, and the 
nonlinear Landau-Zener tunnelling is realized in this system. Finally, some conclusions 
and further discussion are given in Section |VT[ 

II. MODEL 

We consider the Bose-Einstein condensate with total particle number N which is trapped 
in a circularly symmetrical trap V or container, and rotated by an asymmetric perturbation 
V ext . The N-body Hamiltonian in the rotating frame of such a system is given by H = 

N 

J2 [H 0n + V ext (f n ) - u ■ L zn ] + U5(n - fj), where H Qn = p n 2 /2m + V(f n ) is the 

n=l ±<i<j<N 

single body Hamiltonian in the laboratory frame, U = 47rah 2 /m, a is the s-wave scattering 
length and a > for repulsive interactions, to is the angular frequency of the rotating frame. 
The lowest-lying S and P states of the Hamiltonian H are separated by an energy hw . 
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Following m, we adopt the two-mode approximation, i.e. under the condition that the 
effective interacting energy Un is much smaller than Hujq and the energy of the first radial 
excitation is much large than ftwo, where n is the particle number density, we expand the 
field operator in terms of the single particle wave functions of the two lowest-lying states 
which are exact eigenstates of L z , ty(f) = ai<f\{r) + a Lp Q (r). Following the conventional 
procedure and omitting constant terms proportional to N, we get the second quantization 
Hamiltonian 

H = hSco(aQ a — afai) — Vo(ao a\+afao) + ~^(9o a o a t flo^o+^i^^^ a i a i+4^oi a o a o a i" a i)) (1) 

where Su = 0Jq), the coupling energy V = — f ip\{r)V ext l -po{r)dP'r is set real and positive 
by fixing a proper phase difference between ip and <pt and it can be adjusted by changing 
the strength of the perturbation. The g's are determined by the eigenfunctions (fo and (pt 
in the following way g = J ^{rj^dPr, g x = J |^i(f)| 4 rf 3 r, g 01 = J |y2 (r)| 2 |^i(r)| 2 (i 3 r. 



III. MACROSCOPIC QUANTUM COHERENT TUNNELLING 

With the help of the mapping 

-(a^ao - a+ai) -> J z , 
2 ( a o~ a i + a t a o) ~^ Jx, 

^(a£a + at) -> J, (2) 

— * 

our problem is reduced to the one of a particle with large spin J in an arbitrary J-dependent 
potential, i.e. the Hamiltonian of the system can be rewritten as when dropping lower terms 
in self-interactions 

H = h5uJ z - V J X + ^[g (J + J z ? + W ~ Jz? + 4~9oi(J 2 - J 2 Z )] (3) 
First of all, we introduce the SU(2) coherent states [||] 
\Q) = | M ) = e ^+-^-)|-,--> 
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E \l jNT J!!n _„v ( cos f ) ' ""[expHfl sin |] ?+». (4) 



x v (f + n)!(f -n)P 2' - - 2 J 



_> 



Here, 



-itj> 



(0<9<n,0<(j)< 2tt), 

N \ 



(5) 



|f, -f ) is the extremal state | J = f, J z = — f), \n) is the eigenstate of J 2 , 2n represents the 
occupation number difference between the two states and <fr stands for the phase difference 
between them which will be given explicitly in the next section. SU(2) coherent states are 
overcomplete in the Hilbert space and any wave function can be expanded in terms of them, 
for example, a Fock state \n) can be written as 

r-2- 



\n) oc 



i(n+N/2)<t> 



(6) 



In the SU(2)-coherent-state path-integral representation, the Euclidean transition ampli- 
tude from an initial state to a final state is [Bl 



(n f \ e - H{T f~ T ' )/h \ Q t ) = J [dfi(r)] exp 



~\ Se{ 



where 



S } 



in—yl — cost 



(7) 



(8) 



The first term in Eq. (Q) is the Wess-Zumino term, and the semiclassical energy E(9, <$) 
accurate to a constant proportional to iV is given by 

—E{e,4>) = [hdu + -{go - gx)\ cos6» - V sin 9 cos <p + - \g Q1 - -{go + gi)) sin 2 6», (9) 

where 



go = Ug N, 
0i = UfrN, 
9oi = Ug 01 N. 



(10) 



Since sin 6 > (0 < 9 < tt), for Vq > the minima of the semiclassical energy can only 
occur for = 0. 

It turns out that the qualitative behavior of E(9, 0) depends on a relation between Vq and 
the g's. The first regime is characterized by the presence of two minima. For 5uj ^ ^{91—90) 
the minima define a ground state and a metastable state. For 5u = jn(gi — go) the two 
minima are degenerate, and a situation of macroscopic quantum coherence occurs. In the 
other regime E as a function of 9 has one minimum, the position of which depends on 8u. 



The ground state is a coherent superposition of ipo and ipi. In this section, we confine our 
discussion in the first regime. 

For #01 - \(9o + 9i) > and V < g i - \{g + 9i), there always exists 5u = ^(g x - g ), 
such that there exist two degenerate minima in the semiclassical energy. In this case, the 
true ground state of the system will be the superposition of the two macroscopic states due 
to the coherent quantum tunnelling between them which may be observed by interference 
experiments in principle. 

As we have calculated explicitly, such conditions can be satisfied by a large class of 
systems such as Bose gas in an annulus or a thin cylindrical container. Tunnelling between 
the two degenerate macroscopic states leads to energy splitting between them. We evaluate 
this tunnelling splitting by applying the dilute instanton gas approximation which is valid 
for large- N. 

By adding some constants, the semiclassical energy can be rewritten as 



E(8,4>) = Ng(sm9- sinfl ) 2 + Wo sin 0(1 - cos0), 
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where g = |[<?oi — \{% + <7i)], sin# = From 5Se(0,<P) = 0, we get the instanton 

solution corresponding to the transition from 6 = 9q to 9 = n — 9 : 



9 = arccosf— cos#o tanh (u;&t)] 

i cot 2 9osech 2 (uJi ) T) 



arcsm 



2[l + cot 2 o sech 2 (w 6 r)] 1 /2 
where u>b = 2g cos 9 /h. The associated instanton action is 



Sinst — N 

and the tunnelling splitting is 



— cos 9 + - In 



1 + cos 9 
1 — cos 6>o 



\2n J sin# Vl + cos #o/ 



77 cos t 



(12) 



(13) 



(14) 



IV. JOSEPHSON TUNNELLING 



For V > g i — |(<7o + <7i) an d arbitrary 8cu, there always exists only one minimum in 
the semiclassical energy. In this section we evaluate the fluctuations for the relative number 
and phase around the ground state firstly. In the case where Vq > #oi _ \{go + gi) and 



5oj = 7e(<7i — go), the classical ground state of the system is 9 = n/2, <f> = 0. Expanding 
the Euclidean action Se(9, 4>) to the second-order around its classical value 6S c i = and 
performing integrations over 0i(r) = <p(r) — 0(r) and #i(r) = 9{r) — 6{r) respectively [|J, 
we get the effective actions for 9\ and 0i respectively 



0? + — (Vb + 2o)# 2 

8V 1 2 V y; 1 



V N 



dr 
dr . 



(15) 



The motions of 9\ and 
number fluctuation 



l(V + 2g) ri ' 2 
>i are approximately harmonic oscillators and we get the relative 



A(iVo-iVi) 
and the relative phase fluctuation 

A<b = 



\ 



Vn 



V -2g 



(16) 




(17) 

In order to show the relationship between the N-body identical particle state and the 
condensate wave function, we rewrite the SU(2) coherent state \9,<f>) as 

1 



\e. 



—(e-^sm-at + cos 9 -4) N \Q). 
/M 2 1 2 07 1 1 



(18) 



In coordinate representation, it corresponds to 

(9 9 \ 

tp(r u ■ ■ -f N ) = Y[ f e~^sin-v9i(fi) + cos-y? (n) j • (19) 

From the underlying concept of the off-diagonal-long-range-order (ODLRO) of C.N. Yang 
H, we have 



d 3 r 2 ■ ■ ■ d 3 r N ip*(r } r 2} ■ ■ ■f N )ijj(f' ,r 2 , • • -f N ) = — $*(f)$(f '), 



(20) 



and $(r) = \^N[e~ l( ^ sin |</?i + cos |<^o] is called the condensate wave function. 

As we have calculated, for Vq sufficiently larger than 2g we have A(N — N{) ~ v^V <C N 
and A0 <C 1. This means that (p% and (po can be viewed as almost two condensates and 
their phase difference is almost fixed for a large N. Thus, the ground state of the system 
satisfies the condition of the Josephson effect. It's easy to transfer the Hamiltonian (|l|) into 
a Josephson pendulum Hamiltonian 



4q d 2 i = 



COS < 



(21) 



with the Josephson frequency y 8gVo\A — cos 2 6. 

It should be noted that, in the above we consider the case of Su = -^{gi — go) hi which 
the semi-classical ground state gives cos 9 = 0. The discussions above can be also applied 
to other cases in which cos 9 7^ and cos 9 changes gradually as Su — j^(gi — go) increases 
until its absolute value reaches 1. 



V. NONLINEAR LANDAU-ZENER TUNNELLING 

Consider a two-level system, its adiabatic eigenenergies Eo(t) and Ei(t) increases and 
decreases linearly with time t respectively, crossing at a certain value of to- Taking account 
of the coupling energy V between the two levels, the adiabatic eigenvalues of such a system 
Eo(t) and E\(t) with wavefunctions ipo{t) and ipi(t) avoid crossing. We know from the 
adiabatic theorem that if the system with coupling is initially in the state i(jo(t <C to) and 
t changes infinitely slowly from t ^ to to t ^> t , then the system will remain in the state 
ipo(t 3> to)- However, if t changes with a finite velocity, the final state of the system will be 
a linear combination of tpo(t ^> t ) and ipi(t ^> to)- This means when the two eigenenergies 
are getting close to each other, the system will undergo an instanton-like tunnelling, which 
is well known as Landau-Zener tunnelling || and can be realized in a lot of systems in 
molecular and atomic physics 0] . The essential idea of the Landau-Zener tunnelling can be 
generalized to the many-body mean- field problem with non-linear interactions 0, in which 
the tunnelling can take place even in the adiabatic case. 

In the two previous sections we study the case where the angular frequency u is fixed 
while it can be dependent of time in real life. Thus, it is necessary to consider the case 
where u varies with time, accompanied by a non-linear Landau-Zener tunnelling. 

The equations of motion of the operators ao and d can be derived from the Hamilto- 
nian ([]]) 

= ^ U + U (y° a v a ° + 2 9oiafai)]ao - Vfai 
ih-^r = + U(giafat + 2g 01 a^ a )]ai - V a (22) 

Applying the mean-field approximation while dropping the fluctuations of the field operators 

( ai ) = (at)* = x 1 VN 
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(a ) = (a+>* = x VN , (23) 
we get a non-linear Landau- Zener equation 



1 x 1 




U y 





Mu; + cZi + Ci(|xi| 2 — |x | 2 ) — V 

— Vb —Ma; + d 2 - c 2 (|xi| 2 - \x \ 



)l 









(24) 



where 



ci = ^(#0 - 25-01) , 

= ^(^l + 2^0l) , 

c 2 = 2^i- 2 5b0 ' ( 25 ) 

l^ol 2 ; I x i 1 2 are the probabilities in corresponding states and |xo| 2 + \xi\ 2 = 1- Our model is 
a more general one than that of Ref. for the influence of the interaction on these two 
modes is different. If h5u varies linearly with time, i.e. h5u = at, we know from the linear 
Landau-Zener effect that an instanton-like transition will occur at times near to the crossing 



point of the unperturbed energy levels. As the role of the parameters d\ and d 2 in Eq. (|24|) 
is to shift the crossing point of the unperturbed energy levels, we have reason to expect 
the dimensionless parameter G = (go — gi)/ go will only affect the transition time. Another 
characteristic dimensionless parameter of the system K = (gi — 2g i)/ (#0 — 2g i) is given by 
the imbalance of the particle interactions. 

The system of non-linear equations is solved numerically. Fig. |I| shows that the transition 
probability as a function of K for t — > +00 is nearly monotonic increasing. The plots 
with different values of the parameter G are almost identical, which implies the tunnelling 
probabilities are not sensitive to G as we have predicted. As we can see, the system is in 
a coherent state of the two states (po and (pi because of the nonzero transition probabilities 
even in the adiabatic case. This can also explain the hysteresis phenomenon of the vortex 
formation. 



8 



0.84 



0.83 




0.82 



0.81 



0.8 



0.79 - 



0.78 



0.77 



076 



0.2 



0.4 



0.6 



0.8 



1.2 



1.4 



1.6 



FIG. 1: Tunnelling probability as a function of the dimensionless parameter K is plotted for 
different values (0.1, 0.2, 0.3, 0.4, 0.5) of the parameter G. The plots for different values of the 
parameter G are almost identical to one another. 

VI. SUMMARY AND CONCLUSIONS 

In summary, we have studied the macroscopic phenomena in rotating BEC by using the 
SU(2) coherent state path-integral method. Our discussion can also be extended to the 
rotating BEC with attractive interactions and BEC in an optical lattice. In the latter case, 
we can expand the field operator in terms of single-particle ground state wave functions 
localized in each well, then the coupling energy Vq is replaced by the hopping energy, which 
depends on the overlap of the local wave functions. Some detailed calculation can be found 
in an earlier paper |4]]. We can see that when Vq is sufficiently small, the phase fluctuation 
will be very large, ~ 2n. According to Eq. (|6]), the system is in a Fock region. It can be 
easily understood physically in the rotating case that when the coupling energy is absent, 
the tunnelling is forbidden by the angular momentum conservation law and the relative 
number fluctuation is suppressed. In an optical lattice, small Vq means sufficiently high 
barrier which also prohibits the tunnelling. When Vq is comparable to the self-interaction 
energy, the system is in the Josephson region as we have shown in section IV. Thus, contin- 
uously decreasing Vq will lead the system to undergo a phase- diffusion process and there is 
a crossover from the Josephson to Fock region as exhibited experimentally recently in Ref. 
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J8j. The physics of such experiments can be understood clearly and quantitatively in our 
SU(2) coherent state frame. 
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